Electronic structure of d-wave superconducting quantum wires 
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We present analytical and numerical results for the electronic spectra of wires of a d-wave su- 
perconductor on a square lattice. The spectra of Andreev and other quasiparticle states, as well 
as the spatial and particle-hole structures of their wave functions, depend on interference effects 
caused by the presence of the surfaces and are qualitatively different for half-filled wires with even 
or odd number of chains. For half-filled wires with an odd number of chains A'^ at (110) orientation, 
spectra consist of doubly degenerate branches. By contrast, for even A wires, these levels are 
split, and all quasiparticle states, even the ones lying above the maximal gap, have the characteristic 
properties of Andreev bound states. These Andreev states above the gap can be interpreted as a 
consequence of an infinite sequence of Andreev reflections experienced by quasiparticles along their 
trajectories bounded by the surfaces of the wire. Our microscopic results for the local density of 
states display atomic-scale Friedel oscillations due to the presence of the surfaces, which should be 
observable by scanning tunneling microscopy. For narrow wires the self-consistent treatment of the 
order parameter is found to play a crucial role. In particular, we find that for small wire widths 
the finite geometry may drive strong fluctuations or even stablilize exotic quasi- ID pair states with 
spin triplet character. 



I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tivity (HTS), the origin of the pairing phenomenon in 
these materials has been the subject of intense debate, 
and is still not clarified. Part of the unusual nature of 
HTS which has hindered theoretical analysis is the short 
coherence length, which allows short-wavelength fluctua- 
tions of various types of local order to coexist with super- 
conductivity. Most probes of the nature of the supercon- 
ducting state have been restricted, until fairly recently, to 
measurements of bulk properties and tunnelling through 
relatively large areas. Following the pioneering work of 
Hess et al.0, it was realized that scanning timneling mi- 
croscopy (STM) could provide an atomic-scale picture of 
the superconducting state, particularly useful when ap- 
plied to inhomogeneous situations like the vortex lattice. 
Measurements of this type were subsequently performed 
on high-temperature superconductors |3| . In the past few 
years, scanning tunnelling microscopy on the surface of 
HTS have compiled a novel and fascinating picture of 
the local electronic structure of a few of these materials 
H S i, H S i, H 113. In the BazSraCaCuaOg system 
(BSCCO), one dramatic implication of these experiments 
is that even relatively high quality single crystals dis- 
play inhomogeneous electronic structure at the nanoscale 
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In parallel to studies of the HTS materials, point con- 
tact spectroscopy has been used to study the electronic 
structure of ultrasmall conventional superconducting is- 
lands . Among many fascinating consequences of the 
nanoscale geometry are number parity effects, in which 
the qualitative electronic structure depends sensitively on 



whether the number of electrons on the island are even or 
odd. More recently, superconducting wires of widths tens 
of nanometers [TJ have also been fabricated. Although 
they have not yet been studied by STM or similar meth- 
ods, this should be technically feasible. 

Isolated nanoscale grains and wires of HTS material 
have not been fabricated to our knowledge. While this 
may prove technically quite difficult to achieve due to 
the complexity of the crystal structure, there seems to 
be no fundamental obstacle in the long run. This ap- 
plies as well to other superconductors thought to mani- 
fest unconventional superconducting order, where effects 
of finite geometry should be easier to see since coherence 
lengths tend to be larger. 

It is our purpose in this paper to study how d-wave 
(e.g., HTS) and other unconventional symmetry super- 
conductors behave in finite geometry at the atomic scale. 
Ziegler et al.^^ began the study of this problem in the 
case of (100) d-wave quantum wires, pointing out the 
dependence of the Fermi level density of states on the 
parity of the wire width. This is a natural consequence 
of the discretization of the electronic energy levels due to 
the finite wire width in the d-wave state. While in the 
s-wave case all the interesting physics is tuned by the 
level discretization, one expects a priori one fundamen- 
tal difference in the d-wave case: for any geometry with 
surfaces making an arbitrary angle with the crystal axes, 
pair-breaking processes take place on a scale of the co- 
herence length. The most important consequence for the 
electronic structure should be the formation of Andreev 
surface states. 

However, little is known about how Andreev surface 
states behave when the size of the superconductor be- 
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comes comparable to ^o- The zero-energy states form 
on surfaces with orientations different from the antin- 
odal directions of a d-wavc superconductor, due to the 
sign change of the order parameter. In high-temperature 
superconductors, such states manifest themselves as the 
zero-bias conductance peak (ZBCP) in tunneling spec- 
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the anomalous temperature behavior of the 
Josephson critical current [i^, l4ll 142 . \^ and the upturn 
in the temperature dep endence of the magnetic penetra- 
tion depth OliHig (see also review articles |47ll48l|). 
The conventional description of Andreev surface states, 
as well as the Andreev reflection itself, is based on the 
quasiclassical approximation, a powerful tool in the the- 
oretical study of various properties of inhomogeneous su- 
perconducting systems. 

The quasiclassical theory of superconductivity gives a 
so-called coarse-grained description of the phenomena, 
averaged over interatomic distances. This has been used, 
for example, to calculate the local quasiparticle density 
of states (LDOS). However, these coarse-grained aver- 
aged results are not adequate to analyze atomic reso- 
lution measurements by STM and some other contem- 
porary experimental techniques (e.g. atomic force mi- 
croscopy |42|). To obtain this type of information, a fully 
quantum mechanical atomic-scale approach going beyond 
the quasiclassical approximation in describing inhomoge- 
neous states of superconductors is required. We address 
this problem in the present paper using a tight-binding 
BCS-like model of d-wave superconductor on a square 
lattice. 

As a first step towards understanding the effect of con- 
strained geometry, we study the simplest case of d-wave 
superconducting wires consisting of N parallel chains as 
the system size is reduced. The quasiparticle spectrum of 
such systems is described both analytically, with an as- 
sumed spatially homogeneous order parameter, and also 
numerically with a fully self-consistent approach. In the 
limit iV ^ 1 the usual surface Andreev states in lattice 
models 0, HE 03 1 ^^^^ the surface states known in 
continuous models , can be recovered at each surface 
of the wire. However, for sufficiently narrow wires, when 
the transverse wire dimension is the order of the super- 
conducting coherence length, the Andreev states strongly 
interfere and give rise to qualitatively new effects. We 
show below that only those effects which occur for bands 
sufficiently far from half-filling and relatively wide wires 
can be described with the quasiclassical theory of super- 
conductivity. In addition, we demonstrate how and under 
what conditions one can recover earlier quasiclassical re- 
sults for Andreev states in d-wave superconducting films 
[soj from our microscopic approach. 

The microscopic method adopted in this paper to con- 
strain the geometry involves introducing lines of impuri- 
ties of potential strength taken to infinity to bound the 
wire. We show that earlier microscopic results on (110) 
surfaces IH H IM El E3 and (100) wires O can be 



reproduced by this technique, and then extend it to cal- 
culate new results on wires with other orientations. We 
find that the results for electronic spectra are very sensi- 
tive to the number parity of the wire width, and that true 
zero-energy Andreev states can only exist in wires with 
odd numbers of chains. In even wires, the Andreev states 
are split, pushed away from the Fermi energy, and can 
have either surface or standing wave character. Finally, 
we show that for smaller wires self-consistency effects be- 
come important and can even, within mean field theory, 
lead to condensation of a fully gapped spin triplet state 
instead of d-wave order. 

Our results may also have some qualitative relevance 
for the related problem of d-wave superconducting grains, 
where the geometry is constrained in 2 dimensions. In 
addition to artificially fabricated islands, some authors 
have proposed that the BSCCO-2212 samples which dis- 
play nanoscale inhomogencity should be thought of as 
a collection of weakly coupled d-wave grains of roughly 
the size of the superconducting coherence length ^o, or d- 
wave grains coupled to grains of another electronic phase 
0,lalifllS2l- In fact, the structure of a general, possibly 
irregular small grain of d-wave superconductor has not 
been studied to our knowledge. Understanding how the 
local density of states (LDOS) of these wires depend on 
the wire width and the orientation, as well as on the devi- 
ation from half-filling, could provide important intuition 
for the question of the electronic structure of the small 
irregular grains possibly present in BSCCO samples. 

The outline of the paper is as follows. In Sec^I we in- 
troduce the formalism for the problem. In Sec lIIII we 
discuss three special semi-infinite surface orientations: 
(110), (210) and (100). In Scc lIVB2l we study the 
(110) superconducting wires with even and odd width 
and make an comparison with the discrete states in nor- 
mal metal wires fSec lIV 1?T|) to try to identify the nature 
of the true Andreev states. In Scc lIV B'3l we study the 
effects of deviations from the half-filling. In Sec0 the 
results of fully self-consistent calculations are presented. 
In particular, we allow the order parameter to vary spa- 
tially and comment on the differences in our results. In 
the case of narrow wires the self-consistent study shows 
the appearance of some peculiar types of superconducting 
pairing. Finally, in Sec lVIL we present our conclusions. 



II. MODEL DESCRIPTION AND FORMALISM 

The Hamiltonian for a pure singlet superconductor can 
be written as: 

(i,j),<J i,(T 

where we have chosen a nearest neighbor tight-binding 
band for simplicity; /i is the chemical potential. A 
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superconducting pairing is defined for nearest neigh- 
bors Ay = —V{cjiCii) on the bond The pa- 
rameter t is of order 150 mcV for high Tc materials, 
and we consider both the particle-hole symmetric model 
/i = and the more realistic case of finite fi. In 
non-self-consistent calculations the OP has the famil- 
iar fc-space form Ak = Ao[cos(fcaa) — cos(fcba)], where 
Aq = ■^J2±i^ii±ra ~ ^iiiri,) is independent of i, and 
is taken to be 0.2t. The lattice constant is denoted by 
a. The maximum gap is A^ax = 2Ao. We also present 
self-consistent calculations, in which case V is chosen to 
yield this same value of Aq far from wire edges. 

It is possible to constrain the geometry underlying Eq. 
in several different ways. We present results here for 
a method discussed, for example, in Rcfs. lB^ - ls^ in which 
the on-site potentials Ui arc chosen to lie on the bound- 
ary and their value is taken to infinity to cut off electron 
transport through the boundary. This technique has the 
virtue that the strength of the barrier can in principle 
be lowered to allow different degrees of transparency and 
the study of tunneling phenomena. In this work we re- 
strict ourselves to impurity configurations and strengths 
for which the constrained system is completely isolated 
from its environment. This is equivalent to the assump- 
tion of open boundaries, when no hopping and no pairing 
take place outside the region. No disorder is introduced 
in the system and we consider Ui as "impurity" potentials 
only so as to form the surfaces of the superconducting re- 
gion. 

The Hamiltonian of the surface term is taken to be 



clci, 



(2) 



where the set of sites £ is determined exclusively by the 
boundaries of the desired system (see below). The full 
Fourier space Green's function for the system in the pres- 
ence of these impurities is quite generally: 

G(k,k',w) = GW(k,c^)(5k.k' + 

+ GW(k,c^)T(k,k',c^)G(°)(k',c^) , (3) 

where the T-matrix can be found from the following equa- 
tions: 

T(k,k',a;) = [/(k,k') + 

+ C/(k, k" )G(") (k" , oj)f{k' , k', ij).{A) 

k" 

Here G and T take 4x4 matrix form in the four- 
dimensional product space of particle-hole and spin vari- 
ables. If we choose nonmagnetic on-site potentials and 
consider singlet superconductors, the problem reduces to 
2x2 matrices in Nambu space. The Nambu retarded 
propagator for the pure d-wave superconductor is 



where the Tq, are the Pauli matrices and = 
—2t[cos{kaa) + cos(fcba)] — fj,. Calculating the T-matrix 
allows us to obtain the eigenenergies of the system from 
the condition detT~^{u!) = 0. 

The local spin-resolved quasiparticle density of states 
is given as 



(6) 



^tq + CkT3 + Akfi 



(5) 



After integration of the LDOS over energy, we should 
obtain the total number of quasiparticle states per one 
site. Since each site on the lattice possesses 2 states with 
opposite spins, the spin-resolved LDOS normalization is 



dwpi(l)(r,w) = 1. (7) 



III. SURFACE CASE 

Upon a conventional reflection on the (110) surface 
of a d-wave superconductor, the order parameter always 
changes sign as the direction of the quasiparticle momen- 
tum k is varied. This leads to Andreev reflection and, 
eventually, to the formation of the dispersionlcss zero- 
energy surface Andreev bound states. For other surface 
orientations the sign change docs not take place for all 
incoming momentum directions. It is important to no- 
tice that the number of consecutive impurity lines that 
are needed in order to cut the system depends on the ori- 
entation of the surface one wants to consider. For (100) 
and (110) surfaces and nearest neighbor hopping, one line 
of impurities is sufficient to cut communication between 
the two sides. For a (210) surface, the nearest neighbor 
hopping and pairing terms can still connect this particle 
with another across a single impurity line, so for a sim- 
ple tight-binding band, two consecutive lines are needed 
to close the system (Fig. QJ. Alternatively, if one in- 
cludes a next-nearest neighbor hopping t' , even a (110) 
surface is not closed by a single line of infinite impuri- 
ties; the system considered in Ref. [H^ is therefore not 
a closed (impenetrable) surface. For a (hlO) surface in 
general, max(/i, I) lines are needed for a model that in- 
cludes nearest-neighbor terms only, and h -\- I lines are 
needed for a model that includes next-nearest-neighbor 
terms. Clearly, the technique becomes cumbersome for 
arbitrary angles, but one can nevertheless learn a good 
deal by considering special cases. 

In the presence of a surface of arbitrary orientation, 
the simplest way of applying Bloch's theorem to this 
discrete swtem is by using a surface-adapted Brillouin 
zone lSGllsij . We define new coordinates {x, y), rotated 
with respect to the crystal axes (a, &), where x is the 
direction normal to the surface and y is the direction 
along the surface. The system is periodic along the y- 
direction and the crystal momentum component ky of a 
quasiparticle is conserved. Instead of the usual square 
Brillouin zone ka = [— 7r,7r], kf, = [— tt, tt] (for unit lattice 
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(h,l)=(2,l) 




FIG. 1: (210) surface or wire. Two lines of impurities (closed 
squares) are needed to isolate the lattice sites (closed circles) 
with nearest neighbor hopping. 



The Fourier transform with respect to of Eq. |31 is 

G(n, n\ ky, lo) {n - n' ,ky,uj) - (14) 

G^°'> (n, ky,uj) \g^°^ (0, ky, uj)] G^°'> {~n\ ky, cu) 

Due to periodicity of the system along the j/-dircction, 
calculation of the LDOS at site n will simply involve 
a sum of G(n,n,ky,Lu) over all values of ky within the 
surface-adapted Brillouin zone and over two spin direc- 
tions: 

2 f^"^ dk 

p{r, u) = p{n,uj) = --Im / --^Gii{n,n,ky,uj) 

(15) 



constant a = 1) we now use the surface-adapted Bril- 
louin zone given by kx = [—Tr/d, -k / d] and ky = [—ird, nd]. 
Here d = 1/ \/hF+P is the distance between the nearest 
chains (layers) aligned along the surfaces. The momenta 
in the two coordinate systems are simply related through 
rotation of an angle 9 = tan^^ h/l. 

Now we turn to the solution of the equation for the 
T-matrix (Eq. 0)) for the case of one line of impurities. 
We start with the ansatz 

f = rofo + Tifi + r3f3 (8) 

and find, for arbitrary strength of impurity potential, 



TQ{ky,Lu) 
Ti{ky,uj) 
T3{ky,uj) 



Gl^\Q,ky,Lo) 

-Gf\0,ky,u;) 

c-Gl'\Q,ky,Lu) 



(9) 
(10) 
(11) 



where L»i {ky ,uj)^[c- G^"^ (0, ky, uj)]^ - G^^'> (0, ky,uj)^ + 



(0), 



Gf\o,ky,u!)'^, c~ 1/Vq and G^^ ' [x , ky , uj) is the Fourier 
transform with respect to kx of the i-th Nambu compo- 
nent of the bare Green's function G^^^ (k, lo) 



(0), 



Gf\n,ky,uj) 



d 

2^ 



Tv/d 



Gf\kx,ky,u)e'''^"''dkx .(12) 



The site index n corresponds to x-coordinate x = nd. For 
the case of infinitely strong impurity potential considered 
here, c = 0. In this case the expression for T can be 
written in the compact form 



f{ky,Lo) = - GW(0,fey,u;) 



and the poles of the T-matrix correspond to zeros of the 
determinant of the Green 's function G(°) (0, fcj,, w). 



A. (100) surface 

Since (100) surfaces are not pairbreaking in d-wave su- 
perconductors, we do not a priori expect to see interest- 
ing physics arising from Andreev states. On the other 
hand, the mere presence of a surface can induce surface 
Tamm bands j54| . decaying in the bulk on the atomic 
scale and experiencing the Friedel-likc oscillations of the 
LDOS. In our model Tamm states have nothing to do 
with the superconductivity, although the possibility for 
pairing of electrons occupying these surface states is not 
excluded js^ . While it is not our intent to study these in 
detail, we present some results to show that such states 
can be seen by the STM even in situations when sur- 
face Andreev states are absent. FigEl shows on upper 
and lower panels the LDOS for various distances from 
the surface in the normal metal and the superconducting 
states, respectively. 



B. (110) surface 

Results for a (110) surface are expected to be quali- 
tatively different from the ones obtained for a (100) sur- 
face. The bulk order parameter in the coordinate system 
of the crystal axes a and h is Ak = Ao(cos fc^a — cos /c^a), 
but in the coordinates of the surface x (perpendicular 
to surface) and y (parallel to surface) becomes Ak = 
2Aosinfca;dsinfcj,d, with d — In this case, an in- 

cident particle with any non-zero ky experiences a sign 
change in the order parameter as it reflects from the sur- 
face. 

Calculation of the bare Green's function &'^\n, ky,uj) 
from Eq. 1121 gives at ^ = 0, Gi^sin, ky,iij) = for even n 
and Go{n, kyjUj) = for odd n. Explicitly, 



(13) G'^°\n = 2m,ky,u;) = 
iluj 



exp{-i|n|zsgn[tj((j2- A^)]} _ 

: —Tq (16) 

V(w2_A2)(q2„^2) 
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(B/A 



FIG. 2: Local density of states p(x,uj) for a (100) surface. 
Normal state (upper panel) and superconducting state (lower 
panel), p[x,u)) vs. u/Amai for various distances n = x/d 
from surface; Amax = QAt, jj, = 0. 
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FIG. 3: Local density of states p{x,ui) in superconducting 
state versus uj/ A,nax for the case of one (110) surface, Ao = 
0.2t, jj, — 0. Left: chains an odd distance n = i/d from surface 
at x = 0. Right: chains an even distance from a; = 0. 



A^-tj^ < 0,thc branch y/A^ 



for either A = A, g. In Fig. Owe show the LDOS spectra 
on the different layers. Each layer is defined as an array 
of sites parallel to the surface. Its index indicates its po- 
sition; layer n corresponds to sites at a distance nd away 
from the surface. 

The Andreev bound states are manifested as zero- 
energy peaks in the LDOS. For = 0, such states are 
found in all odd layers, and are absent on even ones, as 
seen in Figure 0] This even-odd effect can be easily un- 
derstood from the form of the T-matrix and Green's func- 
tions Eqs. (|13I17|) . For even n, the bare Green's function 



and 

G(°)(n 2m+l,ky,uj) = 
i/S. exp [— i|n|z sgn {(^{q^ — A^))] sgn (n(g^ — A^)) 

~ ^(g2_ A2)(tj2_ A2) 

igexp [— i|?T,|z sgn {Lo[q^ — A^))] sgn((jj) 



Tl 



V(g2-A2)(g2_^2) 



7"3 



(17) 



where 



Z{ky) 
A{ky) 

q{ky) 



tan ^ 
2A, 



9 2 

— 



Cj2_ A2 

kosin(fcy(i) , 
Atcos{kyd) , 



(18) 

(19) 
(20) 



and A{ky) and q{ky) are the maximum gap Ak and 
single-particle spectrum for fixed ky in 110 geometry. 
Note that here the square root function takes positive 
values for positive arguments, i.e. under the conditions 
\/\{ky)\ < \to\ < \q{ky)\ or \qiky)\ < \u;\ < \A{ky)\. FoT 




FIG. 4: Amplitude of Andreev state vs. distance n — x/d 
from 110 surface for Ao = 0.2t, p = . 



G^^^{n,ky,uj) is proportional to lu for low- frequencies. 
Then the T-matrix T ~ l/w for small uj. Even though the 
T-matrix has a pole at = 0, the product of two Green's 
functions in the analog of Eq. ^decreases faster, resulting 
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in zero LDOS at zero frequency. As seen from the right 
panel of Fig. O the low-energy density of states on even 
layers near the (110) surface is substantially less than in 
the bulk, where it varies linearly with sufficiently low en- 
ergy. Furthermore, the amplitude of the gap features in 
the LDOS is noticeably suppressed with decreasing dis- 
tance from the surface. These features of the LDOS can 
be understood based on a simple relation between the 
Green's function G{n' = n,ky,uj) of the half-space and 
the bare Green's function &'^^{n — n = 0,ky,uj), which 
takes place for even n = 2m : 

G{n' = 71, fcj,, cj) = |l — 

- exp [~2i\n\z sgn (cj(q2 - A^))] jc^") (0, ky, uj). (21) 

The factor in the braces in Eq. 1)21(1 controls the deviation 
from the bulk behavior and diminishes with decreasing 
even \n\. At low energies only narrow regions of ky near 
the center and the edge of the BriUouin zone contribute 
to the LDOS. As a result, for even layers the density of 
states goes as |a;|^ at low energies, with the main con- 
tributions arising from ky near the edge of the BriUouin 
zone. 

For odd layers, the ti and T3 components of &^\n = 
2m + l,ky,u!) are the ones that are non-zero, and they 
approach a constant value as lo goes to zero. So the pole 
in the T-matrix generates the peak in the LDOS, asso- 
ciated with the zero-energy Andreev surface states. For 
Andreev states, z is an imaginary quantity. The size of 
the peak decreases as the distance to the surface increases 
due to the e"!"^™^! factor in Ea. (|17|l . For large |ri|, small 
ky dominates the integration over ky in the LDOS and 
we obtain the following zero-energy asymptotic behav- 

2t 

ior of the LDOS: p{lj) = — - — :^S{iLj). The size of the 

ttAow"^ 

zero-energy peak in the LDOS cx n^^. 

C. (210) surface 

It is useful to study a case intermediate between the 
standard (100) and (110) surfaces to see what qualita- 
tively new features arise. From the usual quasiclassical 
viewpoint, the weight of the zero-energy Andreev (210)- 
surface states should be finite, but smaller than for the 
(110) surface because the phase space for which the re- 
flecting quasiparticle experiences a sign change of the or- 
der parameter is reduced. The tight-binding model leads 
to a more complicated dependence of the weight of the 
zero-energy states on the surface orientations relative to 
the crystal axes. Thus, for the (210) surface the model 
shows at half-filling no zero-energy Andreev states at all. 
We associate this discrepancy, in particular, with the dif- 
ference between reflection channels incorporated in the 
two approaches. 

Standard quasiclassical considerations imply that par- 
allel to a smooth surface the momentum component ky 



is conserved in a reflection event, and only conventional 
specular reflection takes place. A tight-binding model 
shows that, generally speaking, this is not the case, since 
the crystal momentum component ky can also change in 
a reflection process by a reciprocal crystal vector along 
the surface. Due to a difference between reciprocal crys- 
tal vectors at the surface and in the bulk, the momentum 
acquired by a quasiparticle in a reflection event can be 
physically distinguished in the bulk from that of specu- 
larly reflecting quasiparticle. Hence, specific crystal pe- 
riodicity along a particularly oriented surface can result 
in additional channels for quasiparticle reflection, if there 
is a reflected state kx{ky), ky on the Fermi surface corre- 
sponding to the ky surface Umklapp process [s^ . 

In this case the Fermi surface, considered as a part of 
the surface adapted BriUouin zone, should exhibit multi- 
ple values of outgoing kp^x for a fixed value of kp,y. This 
situation is realized for the comparatively complicated 
multisheet structure of the Fermi surface in the surface- 
adapted BriUouin zone for the (210) surface, as obtained 
in Fig. 9 of Refl36l A strong dependence of the particular 
shape of the surface- adapted BriUouin zone, as well as the 
Fermi surface in the zone, on surface orientation is an im- 
portant characteristic feature of the tight-binding models 
|36l IsH . Within the quasiclassical theory, the shape of 
the Fermi surface is usually considered as independent 
of surface orientations relative to the crystal axes. We 
note, that the zero-energy surface states, rule, do 
not exist within the quasiclassical approach, if multiple 
channels for refiection of quasiparticles from an impene- 
trable surface are assumed. In this subsection, we now 
study the LDOS for quasiparticle spectra obtained with 
the tight-binding model at half-filling for (210) surface. 

Technically, we now need to solve for the Green's func- 
tion in the presence of two impurity lines (Fig. In 
the general case of two parallel impurity lines, we cut the 
crystal at an orientation given by {hlQ) and introduce an 
impurity potential 

t/(r) = ;7o^<5(r-R,) , (22) 

3 

where Rj are the points on the two impurity lines, at the 
location of the boundaries. The first boundary is defined 
to be located at x = with impurity sites Rj = jd~^y 
and the other is parallel and located at x — {N + l)d 
with Rj = {N + l)dx +{c + jd~^)y, giving a total of TV 
free chains. Here c is a shift along y-axis of sites on the 
{N + l)-th chain relative to the sites on the 0-th chain. 
For the special case of the (210) surface, we need two 
adjacent lines at a; = and x = — d, so formally this 
corresponds to the case of N ~ —2. 

The equation for the T-matrix (Eq. 0)) for this impurity 
potential can be solved by choosing the ansatz 

f{kx,k'^,ky) = io + iie*''^+''^''^ +t2e-^''^+i'''*'-^ 

^l^^^{N+l)d(k^-k'^) (23) 

In the limit of infinitely strong impurity potential this 
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gives 



to 



-G(")(0)- 



t2 



1 - G(o) (O)-iG(o) (- - l)G(o) (O)-iG(o) {N + 1) 

G(")(0)-iG(")(iV + l)G'")(0)'' 

1 - G(")(0)-iG(")(-A^ - l)G(")(0)-iG(o)(iV + 1) 

GW(0)-i 

1 - G(o)(0)-iG(o)(Af + l)G(o)(0)-iG(o)(-A^ - 1) 
G(o)(o)-iG(o)(_7v_ i)G(o)(o)-i 



1 - G(o) (O)-iG(o) (- - l)G(o) (O)-iG(o) (iV + 1) 
The Green's function in this case can then be written as: 

n — n 

( G(«)(0) G(")(-iV-l)yY G^-n') \ 
VG(°)(iV+l) G(o)(0) ) (vG(°)(iV + l-n')y ^ 

where we have not written the dependence on ky and m 
expHcitly. 




FIG. 5: Local density of states p{x,uj) versus Lj/A„iaa: for a 
closed(210) non-transparent surface at half-filling model and 
Ao = 0.2f, /i = 0. Each curve corresponds to a chain located 
at a distance x /d to the surface. 



The equation for the bound state on (210) surface is 



det 



G(")(0) 
G(")(-l) 



G(o)(l) 
G(")(0) 



= 



(25) 



Fig O displays the LDOS spectra calculated on several 
layers for (210) surface. There is no zero-energy peak 
observed on any chain at all, for the reasons explained in 
the beginning of this subsection. The peaks in the LDOS, 
seen close to the (210) surface at energies ±0.5A„iaa;, 
originate from the gap features taken for the momentum 
along the surface normal: /S.{kF.x: kp.y = 0) = 0.5A. The 
peak at w = lS.{kp,x^ 0) arises for the homogeneous model 



of the order parameter, while for the self-consistent spa- 
tially dependent order parameter it lies slightly below 
A(fci?^a;,0). It is associated with the surface Andreev 
states and decays in the depth of the superconductor. 
These peaks have been theoretically found first in Rcf. 
l2^ with a continuous model and then also for the con- 
ductance with a lattice model [l^. We notice that the 
conductance spectrum shown in Fig. 14 of Ref i36[ is in 
agreement with the LDOS on the first chain (at x = d) 
in FigEl The variations of the LDOS from chain to chain, 
which accompany a large-scale behavior, are the Friedel- 
like oscillations. 



IV. HALF-FILLED WIRES 

We now consider wires where the second line of impu- 
rities confines the system to a finite width, i.e. we restrict 
ourselves to the cases with surface normal along the (100) 
and (110) directions. Semiclassically, a quasiparticle will 
go through multiple scatterings, bouncing back and forth 
between the two walls. Hence we expect that the inter- 
play between Andreev refiection, taking place due to sign 
change of the order parameter, and the energy discretiza- 
tion, due to finite wire width, to yield novel features. We 
again model the surfaces by introducing impurities on 
the appropriate sites to completely isolate the wires. The 
corresponding Green's function is given by Eqs.|23J). The 
bound state energies are determined by the equation 



det 



GW(0) 

G(")(A^- 



1) 



g("'(-a^- 

G(")(0) 



1) 



= 



(26) 



A. (100) wires 

Ziegler et al.^ investigated the problem of (100) d- 
wave quantum wires in some detail, and discovered the 
existence of a number-parity effect as a function of the 
width iV of a mcsoscopic d-wave wire: a finite total DOS 
is found at the Fermi level for odd N and zero DOS 
(with a full gap in the excitation spectrum, not a d-wave 
like gap) is found for even N, at least for a simple tight- 
binding band at half-filling. The differences between even 
and odd N were shown to survive for more general bands, 
as well. For completeness, we reproduce, using our ap- 
proach, some of their LDOS spectra in Fig. 

B. (110) wires 

1. Normal metal wires 

Based on intuition from the surface case, we would ex- 
pect that Andreev states play an important role in (110) 
wires, with geometry shown in Figure A crucial ques- 
tion which arises in the following discussion is, how does 
one identify a subgap state of true Andreev character? 
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FIG. 6: Local density of states p{x,uj) vs. lo/ /S.max for a (100) 
wire of width A'^ = 4 (left panels) and N = 5 (right panels), 
using Ao = 0.2t, fi — 0. 




FIG. 7: Geometry of 110 wire (sites are filled circles) of width 
A'' bounded by 2 impurity lines (filled squares) . Wire is infinite 
in both y and —y directions. 



By merely measuring the LDOS with an STM, for exam- 
ple, one may see several peak structures, not all of which 
will be related to Andreev reflections at the surfaces. One 
set of candidate states which needs to be investigated first 
is the set of discrete dispersive (with respect to ky) lev- 
els which arise simply because of the finite wire width. 
These are of course present already in the normal state 
wire. 

The fca^-integrated bare normal state Green's function 
for an infinite lattice above Tc takes the form : 



TT/d 



iexp[z|n| arccos(— ^/^(fcj,))] 



,(27) 



where q{ky) is defined in Ea. H2()(l . The full Green's func- 
tion for a (110) wire may then be obtained by solving the 
T-matrix equation as 



Gii{n, n'; ky,Lu) 



2sm[n,mnz] sin[(n,, 



N - l)z] 



y^q^{ky)-uj^sm[{N + 1)2 



(28) 



where n,n' — l,2,...,iV, z = z{ky,uj) = 
cos~^ {-uj/q{ky)), Tirnin = min(n,7i'), and rimax = 
max(n, n'). 

For every value of ky , there generally exists a series of 
eigenvalues which are the poles of the Green's function. 
They may be obtained by simply solving 



sm[{N + l)z{ky,uj)] = , 



(29) 



which yields: 



uj„{ky) = -q{ky) cos 



7V + r 



0,1,...,7V + 1 . (30) 



This gives N branches of solutions for v ~ 1, ...,7V, dis- 
tributed symmetrically with respect to the Fermi level, 



and two special solutions {v ~ Q and iV + 1) with 
bj = ±q{ky) which we refer to as defining the "effec- 
tive band edge" . One interesting feature is the exis- 
tence in the normal metal wire of dispersionless zero- 
energy quasiparticle states (ZES) which form the branch 
ly = {N+l)/2 when the width N of the wire is odd. Since 
the group velocity vanishes for dispersionless states, they 
are always localized. We note that the case TV = 1 has 
only the obvious dispersionless ZES since no transport is 
allowed with only nearest neighbor hopping. In the gen- 
eral A'^=odd case, it appears that there is always exactly 
one such localized state (doubly degenerate in particle- 
hole space) for any fixed ky. A deviation from half filling 
shifts the zero-energy states to —fi. Thus, for positive fi 
they are the hole states, while for /x < - electron states. 

Normal state Superconducting state 




FIG. 8: Local density of states for a (110) wire of width TV = 
5, with Ao = 0.2t, n — 0. Left: normal state, p{x,uj) vs. 
uj/Amax', right: superconducting state, p{x,uj) vs. uj/l\max- 
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The contribution of each of these states to the LDOS 
can be estimated by examining the residue near the pole, 
where the Green's function can be approximated as: 



G(n, n'] ky,uj) 
with the residue given by: 

Q,y{n, n) 



UJ — ijJu{ky) ' 



(31) 



N+l 



N+1 



N + l 



where n is the index of the layers. Note from H32(l and 
the spectral representation of the Green's function near a 
pole, it is easy to read off the quasiparticle wavefunctions 
as i^yin) = y/2/{N + 1) sm[mTi^/{N + 1)], i.e. just the 
wave functions of a free particle confined to a box of 
width {N + l)d. 




FIG. 9: LDOS for a (110) wire of width N = 4 with Ao = 
0.2£, fi = 0. Left: normal state, p{x,uj) vs. co/Amax', right: 
superconducting state, p{x,uj) vs. uj/A,nax- 



ky) due to the finite width of the wire and transforms 
into conventional continuous quasiparticle spectrum in 
the massive normal metal in the limit oo. The 

momentum resolved LDOS in the normal metal wire with 
discrete dispersive states takes the form 



N 



p{n, ky, ^) = ^ Q„in = n')S{uj - ujy{ky)) 



(32) The integration over ky gives the LDOS: 



Q 












dky 





(34) 



(35) 



where the sum is taken over those ky and which satisfy 
the equation ujy{ky) = uj. 

The position of peaks in LDOS are then determined 
by extrema of the dispersive branches (|30|) . taking place 



at ky — 0: 



-At COS ■ 



N + l' 



,N 



(36) 



The peaks corresponding to the normal metal wire 
states are seen clearly in Figs. |H1 and O at the eigenfre- 
quencies given by Eq. It is easy to check that the 

weights agree with Eq. H32|) . The LDOS for the normal 
metal wires qualitatively differs from the LDOS for bulk 
normal metals with the nearest neighbor hopping on the 
square. The zero-energy peak (the Van Hove singular- 
ity) in the bulk metal is symmetric as a function of the 
energy and its log singularity is much broader than the 
i5-like peaks we find here for wires. 



2. Superconducting wires 



It is easy to check that for any n, n' the residue Q^, van- 
ishes everywhere for the "states" at the effective band- 
width 1/ — or N + I; we therefore do not discuss them 
further, but focus on the N branches with finite residue. 
For these states the residue can vanish locally on chains 
with number n, for which the quantity nv/ (N -I- 1) is an 
integer. The odd TV ZES, for example, has a residue 

2(sin^)2 

Q (N+i) in,n) — — — j — . (33) 

2 iV + 1 

As is seen from Eq.(j3SJ, the probability density of the 
ZES oscillates with a period 2d, taking finite values only 
for odd n and vanishing on all nearest neighbor sites 
( where n is even ). This is also valid for the states with 
energies ±/j, in the case of the deviation from half-filling 
and ensures no transport with nearest neighbor hopping 
for quasiparticles with energies ±pL. 

All these states are not surface bound states, as one can 
easily check that their amplitude docs not decay across 
the wire. The quasiparticle spectrum is discrete (for fixed 



The T-matrix and Green's function equations are nec- 
essarily more complicated in the presence of supercon- 
ductivity, but they are still tractable in the (110) case. 
The bare Green's functions G^°) are given by Eqs. H16|) 
and (|17|l . as before. One must then solve Eq. H2fci|l . 
which applies to any situation which requires two lines 
of impurities, for the eigenenergies lo. Similar to the nor- 
mal metal case, there are special solutions of Eq. (|26|l : 
Lu = ±A{ky), Lu = ±q{ky) for any N and lo = for even 
N. They are poles of the T-matrix, but do not corre- 
spond to the poles of the full Green's function. 

The behavior of subgap surface states on the narrow 
superconducting wire differs qualitatively from the case 
of the superconducting halfspacc due to the interference 
of the wave functions of the states on both surfaces. Since 
the zero-energy peak in the LDOS for the half-filled sur- 
face vanishes on each even chain ( see Figs.lO, (0J), the 
spectrum of Andreev states on the (110) wires becomes 
strongly dependent on the parity of the number N of 
chains in the half-filled wire. This effect is quite pro- 
nounced for wires whose width is the order of or less than 
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the superconducting coherence length. As we demon- 
strate below, some new qualitative features arising in the 
superconducting state of the (110) wires in the quasipar- 
ticle spectrum above A{ky) ( see Eq. (|19ll ') can also be 
strongly dependent on the parity of N. 

Odd N. For a wire with odd N, only the Green's func- 
tions Hlt)|l with even arguments n — 0, ±(A^ + 1) enter 
Ea. (|26|l . Since at small frequencies these Green's func- 
tion are oc w, it is straightforward to show that the only 
subgap state is a dispersionless ZES lo = 0. Dispersive 
modes exist as well and take the form 



N = 11 



N + 1 



A'^{ky) sin^ 



1, 



N + 1 ' 
N - 1 



(37) 



This exactly coincides with the quasiparticle spectrum 
in a bulk two-dimensional superconductor oj'^iky^kx) = 
^"^{ky, kx) -\~ A'^{ky, kx), if the discrete values of momen- 
tum component across the wire k^^u ~ 'nv/{N -f l)d 
are introduced. We note that the dispersionless zero- 
energy Andreev states are the only true subgap states in 
the spectrum, since the energy |w,y(fcj^)| of the dispersive 
states (|37|l lie above the respective value \A{ky,k,j.)\ of 
the order parameter, for any v and ky. 

There are, as in the normal state case, iV — 1 dispersive 
modes in addition to the ZES. They are doubly degener- 
ate due to the particlc-holc symmetry. This simple result 
can be understood as follows. For a fixed ky the problem 
reduces to a one-dimensional problem. The correspond- 
ing two-component Bogoliubov-DeGennes wave function 
takes its values on sites, resulting in 2N degrees of 
freedom in the system. With this point of view it ap- 
pears natural that for fixed ky the total number of levels, 
which are twice degenerate, is A^. The set of levels with 
positive energies for TV = 11 is represented in panel (a) 
of FigHni 

The LDOS for the superconducting wires with N=5 
is shown in the right panel of Fig. |H1 The dispersion- 
less ZES results in a pronounced peak at zero energy on 
odd layers. Furthermore, each extremum of the disper- 
sive mode ijj{ky) results in the peak in the LDOS at the 
energy LOk^ ^^^^ ■ One series of peaks is associated with 
the extrema at ky ~ 0. Since A{ky = 0) = 0, the peaks 
lie at the same positions as in the normal metal wire 
(|^ . Although they are irrelevant to the superconduct- 
ing properties of the wire, some of these peaks can lie 
at finite energies below the maximum of the gap func- 
tion Amax = SAq. For instance, the lowest position at 
finite energies of the peaks of this series is 4i sin -j^S-j- . 
This can be both above or below Amax, depending on 
the ratio Amax/t and the wire width A^. In contrast 
with the normal metal wires, the dispersive quasiparti- 
cle modes Eq. (|37|l in the superconducting wires have ex- 
trema also at the e dge of the surface-adapted Brillouin 
zone ky = ±7r/(2(i) |57| . They are shown in panel (b) of 
FigElfor the wire with A' = 11. This leads to additional 




0.2 0.3 
co/t 

FIG. 10: Quasiparticle spectra and the LDOS for supercon- 
ducting half-filled (110) wire with N — 11. (a) Dispersive 
modes {kyd vs. uj/t) for iV = 11, Aq = 0.2t, = 0. (b) The 
blowup of the dispersive spectra near the edge of the Brillouin 
zone, (c) The LDOS for energies less than Amax'- n—2 (solid 
line), n=l (dashed line), where a broadening 5 = .005t has 
been used. 



series of (A' — 1) quasiparticle peaks in the LDOS for the 
superconducting wires, which lie below Amax- 



Speaks 



^Ajri^jx sin ■ 



N + V 



N-1 

„ = l,...,^—. (38) 



Panel (c) of Fig^l displays the series of peaks in the 
LDOS for TV = 11. 
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The dispersive states forming the nonzero low-energy 
peaks (I38f) in LDOS, are not Andreev states. They 
he below A{ky) for ky near the edges of Brillouin 
zone, but they are situated above the bulk gap func- 
tion A(kx{i'),ky). The wave function for any of these 
states possesses the finite current of the probability den- 
sity, while for Andreev states the total probability current 
vanishes. 

While positions of the peaks are determined by the 
extrema of the dispersive energies, their weights in the 
LDOS are controlled also by quasiparticle wave functions, 
which form the quantity Q{n, n) in Ea. H35|l . For this rea- 
son the weights of the peaks can substantially differ on 
different layers and may vanish on some of them, quite 
analogously to the normal metal wires (see Ea. l|32() 'l. Due 
to an interference from two surfaces, the wave function, 
taken on even layers, turns out to coincide with the re- 
spective wave function in the normal metal case (apart 
from a normalization constant). On the other hand, the 
wave function on odd layers is a superposition of electron- 
like and hole-like Bogoliubov quasiparticles on the wire. 

Even N. In the case of even N half-filled wires, the 
quasiparticle spectra become more complicated. The 
Green's functions (|16|l . (|17|) with even (n = 0) and odd 
{n — zt{N + I)) arguments enter Ea. (|26|l . which can be 
reduced to the following form 



g(fcj,) tan[(A^ + l)z] = aA(fcj,) tan^: 



(39) 



where a ~ ±1, < z < 7r/2. Solutions z^,aiky) of 
Ea. (|39|l are directly associated with lo, in accordance with 
Eq.dlHll: 

^laih) = I'^ih) COS^ Zy^a{ky) + A^{ky) Sin^ Zy^a{ky) . 

(40) 

Here v = 1, ^nd a are the indices of the solution. 
Comparing Eas. (|40|) . (|37|l for the spectra of odd and even 
wires, one can see that z/d plays the role of effective dis- 
crete values of the momentum component kx (at fixed 
ky) across the wire. Ea. H39|) can be transformed to a 
polynomial equation in tan z of the 7V-th degree, if one 
excludes the special solutions lo = ±A(ky) , ±q{ky) men- 
tioned above. Hence, for a fixed a there are exactly N/2 
positive and N/2 negative solutions for oj, describing N 
dispersive branches of the quasiparticle spectrum. Ex- 
plicit analytical form of the spectra can be easily found 
from Eas. H4U|) . 139|) in the particular cases N = 2 

^c.=±\(q{ky) + aA{ky)), (41) 

and iV = 4 (i^ = ±, a = ±) 

= A2(fcy) + ^q{ky)(q[ky) + aA{ky) 

A{ky) ' 



-5a- 



q{ky) 



1 / 5 — 6a 



q{ky) q^{ky) 



(42) 



Eq. (|42() is defined for all ky in the Brillouin zone for which 
^l.aiky) is positive. Eas. (|^ . describe the quasi- 
particle spectra for non-self-consistcnt wires with small 
numbers of chains. As we will show in Sec. O the self- 
consistent treatment of the problem can lead to impor- 
tant modifications of the results, at least if the width 
of the wire is less than or of order the superconducting 
coherence length. 
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FIG. 11: Quasiparticle spectra and the LDOS for supercon- 
ducting half-filled (110) wire with A'' = 10. (a) Dispersive 
modes (kyd vs. uj/t) for N = 10, Ao = 0.2t, = 0. (b) The 
blowup of the spectra near the edge of the Brillouin zone, (c) 
The LDOS for energies less or equal to A,„aa: : n=2 (solid 
line), n=l (dashed line), where a broadening & — .005t has 
been used. 
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Whereas in the normal metal state, v is the only index 
of the solution (for a given fc^), in the superconducting 
even-A^ wire each branch with fixed v splits into two, 
corresponding to two values of the index a, giving a to- 
tal of 2A^ distinct branches. Since the order parameter 
vanishes for fcy = 0, the splitting is absent in this par- 
ticular case. The splitting is associated with the symme- 
try breaking of the eigenstates with respect to the sign 
reversal of ky. In the normal even- A'' wires, as well as 
in the odd N wires (both in the normal metal and in 
the superconducting states), each separate branch of the 
spectra uj^{ky) is an even function of ky. This is not 
the case, however, for the even N superconducting (110) 
wires. As directly seen from Eq.JSHJ and the relations 
q{—ky) = q{ky), A{—ky) ~ —A{ky), the solutions of 
Ea. (|39|l z^^a{ky) with fixed a are not odd or even func- 
tions of ky , due to the nonzero right hand side of Eq. H39() 
in the superconducting state. Evidently, each separate 
branch of the spectra LiJu.a{ky) possesses the same prop- 
erty. The whole spectrum, however, remains symmet- 
ric with respect to ky — s- ~ky, since the branches with 
positive and negative a simply interchange with each 
other under this transformation. Spectra of the half-filled 
iV = 10 superconducting wire with positive energies are 
shown in panel (a) of Fig.^] They agree with the above 
discussion. For branches with higher energies, the split- 
ting is less than for lower branches. For small splittings, 
5uj (X A/{N + 1). 

The splitting manifests itself by slightly shifting the 
positions of peaks in the LDOS, due to respective shifts of 
extrema of dispersive quasiparticle energies. In contrast 
with the double number of extrema, the total number 
of the peaks in the LDOS does not change. Any peak 
in the LDOS originates from the two extrema, situated 
symmetrically with respect to the sign of ky. Each pair 
of extrema belongs to two spectral branches, which differ 
only by their index a = ±1. If we disregard the splitting 
of the quasiparticle spectra, the positions of the peaks 
in the LDOS coming, for example, from the edge of the 
Brillouin zone, are as follows 

peaks = ±A.max con j^^—^, J/=l,...,y. (43) 

The splitting shifts the peak positions towards slightly 
lower energies |w|, as compared with those in Ea. (|43|l . 
The respective extrema of dispersive energies are shifted 
towards lower values of \ky\ from ky = ±7r/(2(i). These 
low-energy peaks in the LDOS are shown in panel (c) of 
Fig. HI] 

The other peaks in the LDOS come, neglecting the 
splitting, from the center of the Brillouin zone. Since the 
order parameter vanishes in the center of the zone, the 
positions of the peaks are still the same as in Ea. H3(i|) 
for the normal metal state. The splitting, taking place 
in the superconducting state, slightly shifts the extrema 
of the spectra from the center of the zone, so that the 
peaks move to a larger energy values as compared 
with those in Ea. H36|) . 



The LDOS for the superconducting wire with A'' = 4 
is shown in the right panel of Fig|51 In addition to the 
peaks at finite energies, discussed above, there is also a 
well pronounced zero-energy peak there (see also panel 
(c) of Fig. II 1() . This peak originates from the extrema 
of two dispersive branches of Andreev states, which take 
place at ztky^od = ±tan~^(2i/Ao), where the relation 
q{ky.o) ~ ^{ky,o) holds. Indeed, as follows from Eas. (|^ . 
(|18|l and, in particular, directly seen from Ea. (|42|l . the 
energy and its derivative equal zero for ky = zLky^. It 
turns out that the multiplicity of zeroes of the lowest 
dispersive branches of states at the extrema ky = ±fcj/,o 
is A^/2, i.e. uj (x \ky — ky^ol^^"^ in the close vicinity of 
kyfi. This follows directly, for instance, from Ea. H50|) 
given below. We notice, that the lowest dispersive curve 
in panel (b) of Fig. manifests very slow change of 
energy when it touches the zero-energy line. This agrees 
with our expectation to (x \ky — fcy.ol^ for the (110) wire 
with = 10. The respective peak in the LDOS of wires 
with large even A^ diverges at a; = 0, but is not a delta- 
like peak as it is for wires with odd A^. 

The identification of the Andreev nature of the quasi- 
particle bound states in confined geometries turns out 
to be a nontrivial problem, at least in the case of half- 
filled (110) wires with even A'^, where there is a symme- 
try breaking with respect to ky — > —ky. As seen from 
Eas. H39|l - H42|) . the asymmetry can be associated with the 
order parameter behavior A(— fcy) — — A(fcy) and with 
the sensitivity of the quasiparticle energies to the 7r-shift 
of the order parameter phase. A dependence of the dis- 
persive energy curve on the order parameter phase is usu- 
ally an intrinsic feature of Andreev bound states only. We 
consider the vanishing of the probability density current 
as an defining property of Andreev bound states. These 
states should be able to carry finite electric current, how- 
ever. In odd A^ wires only zero-energy dispcrsionless 
quasiparticle states satisfy the above requirements. In 
even A^ (110) wires the zero-energy dispcrsionless states 
do not arise at all. The interference of wave functions 
located near two surfaces of the even N wire induces 
dispersive Andreev branches ±a;^(fcj,), which transform 
into the zero-energy states in the limit of infinitely large 
A^. Moreover, all quasiparticle states in the half-filled 
(110) superconducting wire with even A^, even those lying 
above the gap, turn out to satisfy the above mentioned 
conditions, i.e. possess the properties of Andreev bound 
states. As an example, we describe below in detail the 
structure of the lowest dispersive quasiparticle branch, 
which we denote as Andreev branch tOAiky), although its 
energy varies with ky over a large range, including both 
the subgap and the supergap regions (see, for example. 
Fig. nil). 

For those ky where min{|A(A:y)|, (/(fcy)} < |wA(fcy)| < 
max{|A(A:y)|, (/(fcy)}, the wave function of the state can 
be written on odd layers (n = 2m -I- 1) as: 

(;;).Csgn..sin[(A. + l-n)z]^^J (44) 
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and on even layers n = 2m 



(:;(- 1)^/2 sin(nz) 



1 



(45) 



-isgn/cyy 

In the range of ky for which |u;yi(fcj^)| < 
min{|A(A:y)|, ^(/cj,)}, the quantity z, entering 
Ea. (|40|l . becomes imaginary. Under the condition 
|w^(fcy)| < \A{ky)\ < q{ky) thc wave function takes the 
following form on odd layers (n = 2m +1): 

^) =Ci(-l)™sgn^^sinh[(7V+l-n)zi](^^J^^J (46) 

and on even layers n = 2m 

»)^C,(-ir.,i„h(„,0(_.4,J . (47, 

Analogously, if |Li;A(fcj^)| < q{ky) < \A{ky)\, the wave 
function on odd layers {n = 2m + 1) 

, C2SgnujASmh[{N + l-n)zi]( . \ ) (48) 
vj \isgnkyj 

and on even layers n = 2m 



6*2 sinh(n2 



(49) 



1 

-isgnky^ 

Here C, Ci and C2 are normalization constants, zi = 
\lmz{uiA{ky),ky)\, z is defined in Ea. p8|) and taken at 
LU = UJA{ky). 

The condition |u(n, ky)\ = \v{n, ky)\, which is valid for 
all solutions Eqs. H44() - (|49|l . results in zero total proba- 
bility current density, while the electric current does not 
vanish for the given branch. This ensures the Andreev 
character, as defined above, of all the states in even N 
wires, regardless of whether their energies lie above or 
below the gap. 

The quasiparticle states which belong to the same 
dispersive branch ujA{ky) can have their wavefunctions 
both symmetrically decaying in the depth of the wire 
(Eqs. H46|l - H49|) ') or oscillating and forming standing 
waves across the wire (Eqs. (I44II . (|45|) '). In particular, 
the amplitude of the wave functions in Eqs. H4()|) - H49|) 
takes its maximum value on layers n = 1 and n = N , 
and manifests Friedel-like oscillations, as the layer index 
changes from odd to a neighbor even value or vice versa. 
On a larger scale the amplitude decays in the bulk of the 
wire symmetrically with respect to two surfaces. 

It is instructive to follow how the above results trans- 
form into the well known dispersionless zero-energy An- 
dreev surface states in the limit of large N . For suffi- 
ciently large width of the wire compared with the co- 
herence length, and for those ky where the wavefunction 
decays in the depth of the wire, the dispersive energy 
u)A{ky) takes a relatively simple form: 

2A{ky)q{ky) 



X exp 



^\q^iky)-A^ky)\ 

-{N + l)sinh" 



l{\A{ky)\,q{ky)} 



^\q^ky) ~ AHky)\ 



It follows from Ea. ljSOII . that the energy of the Andreev 
states is exponentially small and vanishes in the limit 
of infinitely large A''. With decreasing energy, the range 
of ky where the wave function oscillates (see Eqs. (|44|l . 
(|45|l ). converges to the center and to the edges of the Bril- 
louin zone and finally collapses to the respective points. 
Hence, in the limit of very large N the amplitude of the 
wave function decays inside the wire for practically all 
values of ky. This means that for A^ 00 the dispersive 
branch ujA{ky) transforms into the zero-energy disper- 
sionless surface states situated near the two surfaces. 
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FIG. 12: The weight of the zero-energy peak in the LDOS, 
taken at the first layer, as a function on odd (open circles) 
and even (filled circles) A''. 



Fig. 1 121 shows A-dependence of the weight of the zero- 
energy peak in the LDOS. For odd A^, the weight dimin- 
ishes with increasing A^, while for even N it increases. 
In the limit of large A^ the two curves will converge to 
the weight for the zero-energy surface states when the 
surfaces are infinitely far apart. The odd-even effects 
become negligibly small only for Nd ^ ^o- One could 
try to recover the quasiclassical results for wide wires 
Nd 3> a after averaging over odd and even film widths. 
The boundary conditions for quasiclassical propagators 
are taken somewhere on a distance I from the surface and 
imply some uncertainty about the boundary positions, as 
well as the film (wire) thickness (a <C Z <C ^o)- As seen 
from Fig. ^1 averaging of the weights of the peaks over 
odd and even A^ will strongly reduce their width depen- 
dence. This kind of averaging is much closer (although 
not identical) to the ciuasiclassical results on the width 
dependence of the LDOS for the d-wave superconducting 
film ■ Very recent quasiclassical results [H^, treat- 
ing various wire orientations, demonstrate the appear- 
ance of energy bands of quasiparticle states, in particular, 
for (110) wires. Our microscopic model for high quahty 
half-filled wires with fixed number of chains gives, how- 
ever, only a couple of branches Ea. (|50|l for even A^ wires 
and the dispersionless zero-energy states for odd A^ wires. 
We associate the difference between the microscopic and 
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the quasiclassical results with the particular condition of 
half-filling. The quasiclassical approach implies no sin- 
gular behavior of the LDOS in the normal metal state 
near the Fermi surface, whereas the Van Hove singular- 
ity takes place in the LDOS on the Fermi surface for the 
normal metal state of half-filled infinite square lattice. 
An agreement of our microscopic results with the qua- 
siclassical ones arises in the presence of deviations from 
half-filling (see below Eqs.Illl)),|S3). 

As already mentioned, in the even TV wires the particle- 
hole structure of any quasiparticle states with broken 
ky — > -~ky symmetry satisfies the condition |u(n, = 
\v{n, ky)\. Our picture is that for the states above the gap 
this Andreev particle-hole structure is generated by the 
infinite sequence of "overbarrier" (overgapped) Andreev 
refiections, induced by a sign reversal of the order pa- 
rameter, which the quasiparticles experience along their 
trajectories being bounded inside the wire with impen- 
etrable surfaces. This unconventional feature does not 
take place for the states above the gap in the odd N 
wires, since the two surfaces always result in the standing 
waves across the wire with no important interference ef- 
fects in this case. For negligibly small splitting one should 
consider a superposition of two wave functions, describ- 
ing the two split states. Then the Andreev structure of 
initially nondcgencratc wave functions is lost, since the 
moduli of particle and hole amplitudes can easily differ 
from each other. The Andreev structure of quasiparticle 
wavefunctions can be lost also in the presence of devia- 
tions from the half-filling, if ^ is larger or of the same 
order as the splitting. Since the splitting vanishes for 
ky = 0, one can expect that for sufficiently small ky the 
Andreev structure of the wave functions will be destroyed 
even for small fi. In the next section some other conse- 
quences of deviations from half-filling are considered. 



3. Deviations from half-filling 

The shape of the Fermi surface depends on fi and 
has a strong infiuence on the low-energy quasiparticle 
spectrum. Consider, for example, low-energy quasipar- 
ticle states under the conditions \uj\, |A(A:j^)| <C q{ky). 
This ensures that the quasiparticle energies lie close 
to the Fermi surface and one can take their energies 
in the normal metal state to be in the linear form 
oc vp ■ (k — ki;-). Under this approximation, effects of 
the particle-hole asymmetry are small and one can use 
a quasiclassical approximation, which is valid for quasi- 
particles close to the Fermi surface. Then only the or- 
der parameter on the Fermi surface l^{kp^x,kF,y) en- 
ters the equations. For the wire geometry the momen- 
tum component kp^y is an independent parameter, while 
kp^xikp^y), and A{kp^x{kp^y), kp^y) are actually functions 
of kp^y. For ^ = the Fermi surface for the (110) wire 
is a square with sides parallel to a; or y axes. Hence, 
kp,x — -in: I {2d) actually does not depend on kp,y in 
this case and A(±7r/(2d), fcj^-.j,) = ±/S.{kp,y)^ in accor- 



dance with Ea. H19|l . For finite /i we find /S.{kp^xi kp^y) 

A{kp,y) 



<l{kF,y] 



-yjq^{kp,y) - A'^. 



In the case of superconducting wires the equation 
for quasiparticle subgap energies near the Fermi surface 
\u)\ < \lS.{kp^xTkp^y)\ takes the form 



uj^ sinh^ 



^yA'^{kp,x,kp,y) -g;^ 

\Vx,f{kp,y)\ 



+ sin^ cj){kF.y) = A^{kp^x, kp,y) sin^ ^kp^y) , (51) 



where 4){kp^y) = kp^xd{N -\- 1). The lowest branches of 
quasiparticle spectra, which follow from Eq. (|51|) for N = 
24, 25, 26, are shown in FigHS The solution of Eq. lO 
reduces to a simple form in the limit of large TV: 



uJA{kp,y) = ±2 A{kp^y) sin (j){kF.y) 



exp 



-iN+l)d- 



\MkF,y)\ 

Voojikp,y)\ 



(52) 




FIG. 13: The lowest energy branches for N=24, 25, 26. The 
parameters are t — 2.5, fi = 0.2t = 0.5, Ao = 0.2t = 0.5, 

Amax ~ 2Ao = 1. 



For half-filled wires, when /i = 0, the phase (p{kp^y) 
does not depend on kp,y^ being equal to (fiodd ~ 
for odd N wires and (/)ovon = (m -t- l/2)7r for even N 
wires. This difference between the phases 0odd and (f>ovcn 
plays an important role in forming well pronounced odd- 
even effects in the spectra of wires with odd and even 
numbers of layers. Indeed, for half-filled odd N wires 
Eq. H52|) reduces to the zero-energy dispersionless An- 
dreev states. At the same time, for half-filled even N 
wires Eq. (|52|) describes a dispersive branch of quasipar- 
ticle energies, which coincides with Eq. H5U|) under the 
condition |A(fcj,)| ^ q{ky). 
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N=10 N=ll 




co/t m/t 

FIG. 14: Local density of states for a (110) wire with N=10 
(left column) and N=ll (right column); fi = (upper panels), 
fi — 0.2t (lower panels). On each panel the different curves 
represent various chains. 



For /i 7^ the phase (f>{kF.y) noticeably depends on 

kF,y- 



(j){kF.y) = kF,xd{N + 1) = {N+1) cos^i 




(53) 

Qualitative deviations of low-energy quasiparticle spec- 
tra, shown in Fig. 1131 from the respective spectra of 
half-filled wires (see Figs. E| lll|l are associated with 
the behavior of the phase 0(fcy). The odd-even effect 
in the spectra of (110) wires becomes less pronounced 
in the case of finite as found by Ziegler et al. [l3l| . 
For some values of fi the spectra of odd N and even N 
wires may have no qualitative differences at all. As it 
follows from Eq. (|52|l and, in a more general case, from 
Eq. H51|) . for fJ. ^ additional and strong dispersion of 
the spectra comes from the fcj,-dependencc of the phase 
4'{kF,y)- The larger the N, more oscillations of sm(j){ky) 
take place with varying the kF.y Hence, more cxtrema 
of U!A{kF,y) arise. This results in additional peaks in the 
LDOS, which appear only in the presence of finite fj.. As 
seen from Eq. (|53|l . the phase (j> can considerably vary, 
when the film (wire) thickness varies from Nd to Nd — I 
and a ^ / <C ^0, I ^ Nd. Thus, in averaging the spec- 
trum over the film thickness, a large number of respective 
additional peaks arises, filling the whole low-energy band 
in the (110) wire with edges described by Eq. 15211. This 
is in agreement with the quasiclassical results [E^ [s^ . 

The particle-hole asymmetry is another important fea- 
ture of the spectra. It can be well pronounced for finite ^, 
but lies beyond the quasiclassical approximation. Fig 1 141 
displays the asymmetric LDOS, calculated with Ea. H15(l 
for (110) superconducting wires with iV = 10 and TV = 11 
in the case /i = 0.2t. For comparison, the respective 
LDOS for /i = is also shown. A well-pronounced peak 



at uj = ~fi arises with a deviation from half-filling in the 
LDOS for the odd N wires. We remind the reader that 
in the half-filled normal state (110) wires with odd N the 
dispersionless zero-energy quasiparticle states have been 
found in Sec lIVBll (see Eq.^ with v = {N + l)/2). 
The wave function of these states, as well as the residue 
of the pole-like term in the Green's function Ea. (|33|l . is 
a standing wave across the wire, taking zero values on 
alternating sites. In the superconducting state the zero- 
energy standing wave disappears and the dispersionless 
zero-energy Andreev surface states arise. Their weight 
exponentially decays in the bulk of wide odd N wires. In 
the presence of a deviation from half-filling the energy of 
the quasiparticle states in the normal state odd N wires 
shifts to — /i. These dispersionless states with finite en- 
ergy keep the character of standing waves. Further, in the 
superconducting wires with finite fi the low-energy states 
become dispersive both for odd and even N wires. As 
seen from Fig. (|13|l for the wire with = 25, the branch 
with lowest energy has in this case two extrema. The 
maximal value of |a;| for the states forming this branch, 
lies at ky = and contributes to the peak at w = — /i 
associated with the hole contribution, if /z > 0. Since the 
order parameter vanishes at ky = 0, the respective quasi- 
particle wave function is a standing wave and the peak 
position coincides with that in the normal metal state 
of the wire. The minimal value of |a;| is zero and con- 
tributes to the zero-energy peak. The zero-energy peak 
is associated with comparatively large value of ky, com- 
parable with the size of the Brillouin zone, and the order 
parameter at this value of ky is of order Amax- Thus, 
the zero-energy peak is associated with surface Andreev 
states, whose quasiparticle wave function decays in the 
bulk of wide odd N wires. For narrow wires, the self- 
consistency condition becomes important. As shown in 
the next section, at finite fj. the self-consistency condition 
can lead to more important consequences as compared 
with the case = 0. 



V. SELF-CONSISTENT TREATMENT OF 
(IIO)-WIRES 

In previous sections, the order parameter was assumed 
constant over the whole width of the wire in order to 
allow for analytical solutions. Even for a single (110)- 
surface, however, a self-consistent treatment of the order 
parameter gives rise to interesting effects: the d^2_y2- 
wave order parameter is strongly suppressed near the 
surface and a complex is-admixture (or some other 
time-reversal symmetry-breaking state) is possible, which 
leads to a splitting of the zcro-cncrgy Andreev bound 
state mils,, 22, 2 J, HI 111 .39,, 59, 60, il H H |6|. 
Even larger effects are therefore to be expected for the 
wire fimited by two (llO)-surfaces. Indeed, our self- 
consistent evaluation indicates that for very narrow wires 
a quasi one-dimensional triplet superconducting state can 
replace the conventional d^^-y^ + is-state. For finite 
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chemical potential fi the normal metal state can become 
energetically favorable as the ground state for narrow 
wires, while superconductivity recovers with increasing 
wire width. Under special conditions, even a mixture of 
singlet- and triplet-pairing can occur. It is important to 
recall at this point that one expects mean field theory to 
break down as the ID limit is approached even at T = 0. 
Thus the predictions for various kinds of superconduct- 
ing order mentioned below are to be treated with some 
skepticism as regards quantitative predictions. Neverthe- 
less we view our results as presenting intriguing evidence 
that when surface energies begin to become compara- 
ble to the energy differences between c?-wave and other 
bulk pair channels, strong fluctuations with symmetries 
optimal for quasi-lD system, including spin triplet pair 
fluctuations, will result. 



A. Order parameter 

Self-consistent solutions to Hamiltonian were ob- 
tained by solving the Bogoliubov-de Gennes equations 
for (llO)-wires, 



^-ky -t.ky J \VkyXj \Vky\, 



where Ek^x with A = 1 . . . are the eigenvalues of 
the Bogoliubov-dc Gcnncs equations and its eigenvectors 
Ukyxin) and Vkyx{n) are the coefficients of the Bogoliubov 
transformation: 

Ckyu] = ^{lkyX]Ukyx{n)~il^y^^vl^xi.n)^ 

X 

c-fe,ni = ^{lky\]Vkyx{n)+^l^^^Uk^x{n)] (55) 

A 

Furthermore, S^k and are matrices in the N layers 
of the (llO)-wire. i.e. the layers in x-direction, and arc 
given by (for notation see Fig. [7|): 



(Cfe„ 



cos fcj, ((5n'n+l + <^n'n-l) ^ l^'^n 



k) , = (A+ ,e-''=« -I- A" ,e'^y)5n'n+i 

/ 7in' \ nn' ' nn' / n-t-i 



+ (Awe-'^+A^.e^''") (56) 

The gap values are determined by the following self- 
consistency equations: 



nn-\-l 



-^^E«^''''<^-^."+UCfc,nT) (57) 

ky 



ky X 



To simplify the numerical evaluation of the 
Bogoliubov-de Gennes equations we consider iso- 
lated wires here. For a Hamiltonian on a discrete lattice 
like ^ an isolated wire is equivalent to a wire limited 
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FIG. 15: Narrow wires with a nearest neighbor interaction 
strength of l/=1.1575t, which gives rise to a gap value of 
Ao=0.2f for the bulk system at /i=0. Lower panels: phase 
diagram of the wires with even width (c) and with odd width 
(d). White space denotes A = 0, i.e., the normal state. Upper 
panels: corresponding amplitudes of the d^2_y2 and is-order 
parameters displaying the value in the center of the wire for 
the d^2_j,2-wave case and the value at the edge of the wire 
in the s-wave case. These are the positions where the largest 
values of the respective order parameters are expected in the 
usual d^2_y2 + is-state. 



by lines of unitary impurities, which is the boundary 
condition used for the analytical calculations in the 
previous sections. 

For narrow wires we find a variety of different phases. 
The resulting phase diagram for a nearest neighbor in- 
teraction strength of V^=:1.1575t, which gives rise to a 
gap of Ao=0.2t in a bulk system at /i=0, is displayed in 
Fig. 1151 For wires with widths up to A^=9 we find a new 
phase over a wide range of chemical potentials charac- 
terized by Aij = — Aji (stars in Fig. [TCT c). (d)), which 
is a signature of triplet-pairing with Sz = 0, the only 
triplet component compatible with Hamiltonian For 
fi = the amplitude of A^- oscillates across the wire 
between zero and its maximum value (see Fig. I16f a)). 
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FIG. 16: Quasi one-dimensional triplet superconducting state 
for A'^=8 and ^=0. (a) Bond gap values A^^^j^ (for definition 
see Eq. I58II ) across the wire, (b) Illustration of the one- 
dimensional structure of the superconducting correlations for 
an A''=6-(110) wire at ^ = 0. (c) Density of states starting 
from the outermost layer of the wire up to the middle layer 
of the wire, where for illustrational purposes each layer has 
been shifted by an additional offset of 0.3t. 



indicating a onc-dimcnsional nature of these new triplet- 
superconducting correlations (sec Fig. I16f b)). Although 
the oscillating behavior of remains for finite n, its am- 
plitude no longer vanishes exactly on alternating layers. 
Thus the strict one-dimensionality of the superconduct- 
ing correlations seems to be a feature peculiar to fi ~ 0. 
This new triplet superconducting phase will be discussed 
in more detail below. At first we will focus on the singlet 
superconducting phase with possible d^2_y2 - and s-wave 
order parameters. 

Although we do not understand all details of the vari- 
ability of the phase diagram, a few general trends seem 
clear. The larger the width of the wire the more of the 
usual d,^2_y2 -l-is-phase (filled circles in Fig. lT^ c'l.fd)) is 
recovered. The corresponding amplitudes of the dx2-y2- 
and is-components of the order parameter are displayed 
in the upper panels of Fig. 1151 For narrow wires the 
amplitude of the d3.2_j,2-wave order parameter is finite 
only for small chemical potentials ^ and a pure s-wave 
phase is favorable for large fi. These two phases are sep- 
arated by a normal state region. For larger wire widths 
the range of the (ia;2_,y2-wave phase and the amplitude of 
the d^2_j,2-wave order parameter increase. The ampli- 
tude of the s-wave order parameter, on the other hand, 
decreases and the s-wave phase moves towards smaller 
chemical potentials until it merges with the dx2-y2~waNe 
phase, thereby giving rise to a finite is-admixture near 
the edges of the wire, as expected in analogy to a single 
(llO)-surface. 

The dependence of the amplitude of the dx2 _y2-waNe 
order parameter on the width N of the wire is shown in 
Fig. El The upper panels refer to an interaction strength 
of V^=1.1575i, which gives rise to a gap of Ao=0.2i for 
the bulk system at /z=0. For ^=0 we observe an even- 
odd oscillation in the amplitude of the order parame- 
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FIG. 17: Evolution of the d^2_y2-wave order parameter with 
increasing wire width for two different interaction strengths 
V^=1.1575t (upper panels) and V^=1.7682t (lower panels), 
which give rise to a gap value of Ao=0.2t and Ao=0.4t for the 
bulk system at ^ = 0, respectively. The first column displays 
the amplitude of the c?^2_j^2-order parameter as a function 
of the wire width for two different chemical potentials fi/t=0 
and fi/t = 0.4. The other two columns show the variation of 
the d^2_y2- and is-component of the order parameter across 
the wire for fi/t=0 (center column) and fi/t=0.4 (right col- 
umn). For comparison the reduced magnitude of the bulk 
order parameter at /i/t=0.4 is displayed with dashed lines in 
the right panels. 



ter (sec Fig. I17r a)). which disappears for larger chem- 
ical potentials. The effect of finite ^ on the suppres- 
sion of the dj.2_j,2-wave order parameter is considerably 
stronger in narrow wires than in the bulk system. For 
/i/t=0.4 the d^2_j,2-wave order parameter vanishes for 
< 10 (squares in Fig. I17f a)). although the bulk gap is 
only reduced by a factor of approx. 0.8 (see Fig. I17r c)). 
The dramatic suppression of the d^2 _y2--wave order pa- 
rameter for finite n is reduced upon consideration of a 
larger nearest neighbor interaction strengths, as can be 
seen in the lower panels of Fig. El Foi a larger inter- 
action strength of y=1.7682t, which gives rise to a bulk 
gap value of Ao=0.4t at ^ = 0, a pronounced even-odd 
effect remains at ^/i=0.4 (see Fig. Eld)) and the ampli- 
tude in the center of an A^=12-wire is already very close 
to the bulk value (see Fig. ll7r f)) contrary to the smaller 
interaction strength, where the amplitude in the center 
of an A^=12-wire is still suppressed by a factor of more 
than three with respect to the bulk value (see Fig. lTTT c)). 

The emergence of a new triplet-superconducting phase 
for very narrow wires can be most easily understood by 
considering the smallest wire, i.e., the A^=2 wire. Al- 
though the effect of fluctuations will be larger for smaller 
wires, we focus only on possible solutions of the BCS 
mean-field Hamiltonian in the present paper. Inspec- 
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tion of the phase diagram of Fig. shows that the d- 
wave order parameter vanishes for the A'^=2-"wire. To in- 
vestigate alternative ways in which the N^2 wire could 
lower its ground state energy, we map it to a ID-chain 
(see Fig. I18|l . Although the gap values and Aji could 
in principle differ in amplitude and by a phase factor 0, 
we restrict our treatment to (j) — and (j) = n. Note that 
= corresponds to a singlet pairing state, whereas 
= TT is a triplet pairing state with Sz = 0^^- To inves- 
tigate the possibility of triplet superconductivity in these 
systems in more generality, one should retain pairing cor- 
relations with S'z = ±1 in the mean- field Hamiltonian as 
well. For now, however, we are satisfied with the obser- 
vation that even for the simple nearest-neighbor pairing 
Hamiltonian |^ a triplet order parameter can be favored 
over a singlet order parameter in narrow geometries. 



FIG. 18: Mapping of the N—2 wire with (110) orientation to 
the ID-chain with lattice constant d. Note that A2 — — Ai 
corresponds to a triplet pairing state with Sz=0- 

After Fourier transformation Hamiltonian for the 
ID-chain reads 



Hmf = - y^(2t COS fed + ^i)c\^Cko 



ka 



^{(Aie*^'^ 



^26 



(58) 
h.c.} , 



where A^j+i = -F(c,+ixCj| ) = Ai and A^+h = 
— y(ci^Ci-)-ii) = A2. It can be easily diagonalized us- 
ing the Bogoliubov-transformation to give a quasiparticle 
dispersion of: 



Ek 



with Cfc 



qik) 



(59) 



fi = 2t cos kd + fj. 



Al = Al + Al + 2Ai A2 cos 2kd , 

where the gap values are to be determined from the fol- 
lowing self-consistency equations: 



Ai = -V 

A2 = -V 



1 

And 
1 

And 



dk- 



A2 cos 2kd 



Aicos2fcrf + A2 
dk . (60) 

TTd ^k 



Four different phases emerge from this model by varia- 
tion of the chemical potential /i and the nearest neighbor 
interaction strength V (see Fig. I19|) . Below a critical in- 
teraction strength the order parameter vanishes and the 
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FIG. 19: Phase diagram for ID-chain as a function of chemical 
potential /i/t and nearest neighbor interaction V/t. 



normal state is the ground state of the ID chain. For in- 
termediate interaction strengths and small chemical po- 
tentials we find Ai = — A2, i.e., a pure triplet supercon- 
ducting state, whereas for large chemical potentials the 
ground state is characterized by Ai = A2, i.e., a sin- 
glet extended s-wave state. In the singlet state the gap 
function is Afe ~ cos kd, i.e., it has nodes at kd ~ ztn/2 
whereas the triplet state gap function is A^ ~ sin fed 
and therefore is maximum at kd = ±7r/2. This explains 
why the triplet state is favored over the singlet state for 
/i = 0, where the Fermi-surface of the nearest neigh- 
bor tight-binding model is at kd = ±n/2. For larger 
chemical potentials the situation is reversed and the sin- 
glet extended s-wave state becomes more favorable than 
the triplet state. For large interaction strengths, we find 
different magnitudes for Ai and A2, which corresponds 
to a mixture of singlet and triplet pairing. In the limit 
V ^ 00, either of the gap values Ai and A2 approaches 
zero whereas the other goes to infinity, corresponding to 
an admixture of triplet and singlet order parameters with 
equal amplitudes. 

With the phase diagram of the ID-chain in mind we 
arc now able to better understand the phase diagram of 
narrow wires as displayed in Fig. 1151 For small ^ the 
new phase with Aij = —Aji simply arises from the for- 
mation of quasi one-dimensional triplet-superconducting 
correlations along A^=2 wires (see Fig. I16r a).fb)). This 
also explains why this new phase is more favorable for 
wires with even width than with odd width (compare 
left and right panels in Fig. I15|l . as only the former can 
be divided evenly into N~2 wires. Note that also for 
the narrow wires, a critical coupling strength of similar 
magnitude as in the ID-chain is necessary to induce the 
triplet superconducting state, which is, however, smaller 
than the interaction strength considered in this paper. At 
larger chemical potentials we find an extended s-wave or- 
der parameter analogously to the one-dimensional chain. 
More difficult to reconcile with the phase diagram of 
the one-dimensional chain, however, are the seemingly 
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N=10 N=ll 




FIG. 20: Density of states for the A''=10-wire (left panels) 
and the A''=ll-wire (right panels). Upper panels show the 
non-selfconsistent results assuming a constant d^2_y2 -order 
parameter of Ao=0.2t, whereas the lower panels display the 
density of states resulting from the self-consistently deter- 
mined values of the d-wave order parameter, whose variation 
throughout the wire is depicted in the insets. All panels show 
the density of states starting from the outermost layer of the 
wire up to the middle layer of the wire, where for each layer 
towards the center of the wire an additional offset of 0.3t has 
been used. 

arbitrarily distributed mixtures of singlet- and triplet- 
superconducting order parameters (crosses in Fig. I15|l . 
The normal state regions, which occur for narrow wires 
and intermediate chemical potentials in Fig. 1151 reflect 
the fact that superconductivity becomes less favorable in 
finite geometries and a critical coupling strength is nec- 
essary to induce it [s^l • 

B. Density of states 

What does the phase diagram in Fig. El imply about 
the possible existence of Andreev bound states in nar- 
row wires? In Fig. llGf c) the local density of states and 
the variation of the order parameter throughout the wire 
is depicted for the quasi one-dimensional triplet super- 
conducting state using the A^=8-wire at ^=0 as an ex- 
ample. Obviously, the density of states in this state is 
fully gapped, in analogy to the one-dimensional chain, 
and there are no Andreev bound states. 

Close to half-filling, for wider wire widths, the results 
are more conventional and the effects of self-consistency 
much simpler. In Fig. 1201 the local density of states is 
displayed for the A^=10- and the 7V=ll-wire at n = 0, 
which are characterized by a pure dx2_y2-OTdeT parame- 
ter. Here, the main difference between the self-consistent 
(lower panels of Fig. I20|l and the non-self consistent re- 



sults (upper panels of Fig. I20() is the suppression of the 
magnitude of the d-wave order parameter especially to- 
wards the edges of the wire (see also insets of Fig. I2U|I . 
Whereas the weight of the zero-energy peak is reduced 
in the even-width (A^=10) wire the weight of the zero- 
energy state in the odd- width wire (A^=:ll) is hardly af- 
fected. 



VI. CONCLUSION 

Motivated by recent scanning tunneling experiments 
on the Ba2Sr2CaCu208 systems which reveal inhomoge- 
neous electronic structure on the nanoscale, we have ana- 
lyzed in detail the electronic structure of d-wave quantum 
wires. These wires exemplify the effects of a constrained 
geometry while still being simple enough to allow for an- 
alytical solutions. To impose a restricted geometry we 
use lines of impurities with infinite scattering strength, a 
method which allows to cut arbitrarily shaped objects out 
of the two-dimensional plane. In principle, it is straight- 
forward to extend this method to investigate the effects 
of tunneling between neighboring grains by reducing the 
scattering strength of the impurities and thus lowering 
the potential barrier between neighboring grains. 

New and interesting physics arises from interference ef- 
fects between the two surfaces of the wire when its width 
is of the order of the superconducting coherence length. 
Contrary to s-wave superconductors in finite geometries, 
the surface pair breaking plays an important role. In 
this respect, the existence and nature of Andreev bound 
states in constricted geometry is of particular interest. In 
order to single out new effects peculiar to quantum wires 
and arising from the interference of the two surfaces, we 
have first addressed the case of a single surface in the 
first part of the paper, concentrating on surfaces with 
(100)-, (210)- and (llO)-orientations in a nearest neigh- 
bor tight-binding model at half-filling. Andreev bound 
states can form on surfaces with orientations deviating 
from the (100) direction due to the sign change of the 
d-wave order parameter. However, in the presence of 
several channels for reflection of quasiparticles from the 
surface, the zero-energy Andreev states may not exist, 
as is the case for the (210) surface of the square lattice. 
Here our results are in qualitative agreement with earlier 
work based either on the quasiclassical approximation 
or Bogoliubov-dc Genncs equations. Only for the (110)- 
surface, which involves the strongest pair-breaking, do 
we find a zero-energy Andreev bound state, whose am- 
plitude decreases with the square of the inverse distance 
from the surface and vanishes on even layers. 

In the main part of this work, we focussed on quan- 
tum wires with (llO)-orientation, which display a pro- 
nounced "width parity" effect. The special case of elec- 
trons hopping on a square lattice with half-filled tight- 
binding band was treated most extensively. For wires of 
this type with an odd number of chains N only one sub- 
gap state, a dispersionless zero-energy Andrecv-bound 
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state, exists. As in the single surface case, the amplitude 
of the ZES vanishes on even layers and decays towards 
the center of the wire. In addition to the zero-energy 
state (ZES), there are N-1 dispersive modes, which are 
doubly degenerate due to particle-hole symmetry and are 
not Andreev states. 

In the case of even half-filled wires, a splitting of 
the branches occurs which is associated with a symme- 
try breaking ky ^ky, and a total of 2N dispersive 
modes exist. Although there is no dispersionless ZES, all 
quasiparticle states in the even width wire are of Andreev 
character in the sense that the current of the probability 
density vanishes due to opposite contributions from par- 
ticle and hole excitations. These quasiparticles occupy 
either conventional Andreev-type surface states, or a new 
type of Andreev standing wave, according to their mo- 
mentum ky parallel to the wire. In the 2D limit N oo, 
the lowest energy dispersive state was shown to trans- 
form into the usual zero-energy Andreev bound state at 
the impenetrable surface, while the standing wave states 
evolve into either the continuous spectrum or the surface 
states. With increasing deviation from half-filling, odd- 
even effects in the wires become less pronounced. The 
evolution of the Fermi surface shape with these devia- 
tions can result in additional extrema in the quasiparticle 
dispersive modes and, hence, new peaks in the LDOS. 
Large-scale faceting of the surfaces with characteristic 
scales larger than the wire thickness will not influence 
our results significantly. However, small-scale inhomo- 
geneities like point-like defects and impurities can sub- 
stantially change the effects of interference induced by 
wire surfaces, even if the phase breaking length is large. 

The small coherence length of high-temperature super- 
conductors leads to further restrictions on the applicabil- 
ity of the quasiclassical results to superconducting wires 
or films, whose width is less than or comparable to the 
coherence length. For narrow wires of (llO)-orientation, 
the self-consistent treatment of the order-parameter is 
found to have a large effect. For wires with widths less 
than the superconducting coherence length (up to A^=9 
for the particular parameters used in our numerical calcu- 
lations), especially for the even width wires, a new phase 
characterized by quasi one-dimensional triplet pairing is 
found in the mean-field phase diagram. This new phase 
is fully gapped and characterized by the absence of An- 
dreev bound states. The larger the width of the strip, 
the more the d + is-state, which is expected for a single 
surface with (llO)-orientation, dominates the phase di- 



agram near half filling. With respect to the density of 
states at half-filling, where the order parameter has only 
a d-wave component, the main effect of the self-consistent 
treatment is the suppression of the magnitude of the d- 
wave order parameter, especially near the surfaces of the 
wire. 

It is interesting to end this discussion with some spec- 
ulations on the role of bound quasiparticle states and 
edge effects of this type on the spectra of weakly cou- 
pled superconducting grains as apparently observed in 
STM experiments. Such irregular grains should con- 
tain nanoscale "facets" at all possible angles, so presum- 
ably the most general situation with sizeable particle-hole 
asymmetry and mixture of even- and odd-iV boundary 
conditions will apply. If we first assume that the pair 
interaction and grain size are such that one may ignore 
the triplet states found in the self-consistent treatment, 
we expect the spectra of weakly coupled grains to be 
dominated by the zero-dimensional analogs of the dis- 
persive subgap states discussed here, i.e. there should 
be a wide distribution of boimd state energies depending 
on local geometry of the grain, and visible in the LDOS 
as measured by STM. In this sense we question whether 
the weakly-coupled grain picture is in fact applicable to 
the experiments in question, which appear to see a very 
homogeneous spectrum at low energies in the supercon- 
ducting state. 

A second remark is based on our observation that 
in nanoscale confined geometry, spin triplet fluctuations 
may become more favorable. Such time-reversal symme- 
try breaking fluctuations will clearly lead to local spon- 
taneous currents, an issue which has recently been raised 
again in angle-resolved photoemission studies [g^. Fu- 
ture studies of small grains are planned to address these 
issues. 
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